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Abstract 

The first measurement of diffractive processes in deep inelastic scattering at 
the HERA collider is analysed in terms of a "soft" pomeron exchange model. The 
partonic structure of the pomeron which emerges in this picture is determined using 
the QCD parton model. The important role of the gluonic component in the partonic 
interpretation of the data is particularly emphasized. 
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1 Introduction 



The recent measurements performed by the HI and ZEUS collaborations at the DESY 
electron-proton collider HERA Jl|-0 have shown that there is a significant diffractive 
contribution in deep inelastic scattering. To be precise, diffractive deep inelastic scattering 
is the process 

e + p^e' + X + p' (1) 

where there is a large rapidity gap between the diffractively produced system X and the 
recoil proton p' (or excited proton state). In principle the scattered electron e' and the 
scattered proton p' can be identified and their momenta measured. In this case the cross 
section for the diffractive process (P is characterized in the leading twist approximation by 
two structure functions dF® jdxjpdt and dF® / dxpdt in analogy with the decomposition 
of the unpolarised inclusive ep cross section: 



da D 47r 2 a 2 



dxdQ 2 dxpdt xQ 4 



y 2 \ dF 2 D (x,Q 2 ,x ri t) y 2 dF£(x 1 Q 2 ,x 1Pl t) 



(2) 



2 J dxjpdt 2 dxjpdt 

where the kinematical variables are defined as follows: 

n 2 2 Q 2 PI 
Q = -q x = — y = 

2pq p e p 

t = (p- P 'T P = 9 ( Q2 ** = \ ( 3 ) 

2q(p - p') (3 

and q = p e — p' e i where p e ,p' e ,P and p' are the four momenta of the initial electron, final 
electron, initial proton and the recoil proton respectively (see Fig. 1). 

No accurate determination of t variable has so far been possible for the measurement 
of the diffractive structure function at HERA. Therefore the structure function integrated 
over t 

F^(x, Q\ ,„) = £_ it ^fj"^ (4) 

was evaluated from the cross section (g). The longitudinal component was neglected in the 
original procedure of @, |j] because of anticipated smallness of the longitudinal structure 
function and the measured y values, y < 0.5. 

The main observation is that the structure function F®^ factorises 

F^\x ) Q\x ]P ) = x- p n F((3,Q 2 ) , (5) 

where (3 = x/xp. The exponent n is found to be independent of (3 and Q 2 , and equal to 

n = 1.19 ± 0.06(stat.) ± 0.07(syst.) (6) 
n = 1.30 ±0.08{stat.)±° o il(syst.) (7) 

for the HI and ZEUS measurements respectively 0, M. 



1 



Such a universal dependence is expected in models of diffractive interactions ([!]) based 
on the concept of an exchange of a pomeron - an object which carries vacuum quantum 
numbers |5|, |^j . In that case a virtual photon, emitted by the incoming electron, scatters 
off a colourless object (pomeron), which carries a small fraction xjp < 0.1 of the incoming 
proton's momentum. This kinematical condition leads to a rapidity gap between the 
outgoing proton p', which stays almost intact, and the diffractive system X, produced as 
a result of the virtual photon-pomeron deep inelastic scattering. The function F(f3, Q 2 ) 
in factorization formula (|5p is proportional to a pomeron structure function defined in 
analogy with the proton structure function F 2 (x,Q 2 ). The (3 variable plays a role of 
the Bjorken variable x and, if a simple partonic interpretation of the pomeron structure 
function is valid, gives the fraction of the pomeron's momentum carried by a partonic 
constituent of the pomeron. In a more refined approach the QCD-improved parton model 
can be applied and the partonic content of the pomeron can be studied just as in the 
proton case. 

In essence this is the Ingelman-Schlein ([[5]) model. In this paper we show that the 
diffractive deep inelastic scattering data from HERA are described by such a model. In 
this way we investigate the partonic structure of the pomeron. 



2 The Ingelman-Schlein model 

In this approach the diffractive structure function dF 2 jdxjpdt is written in the factorisable 
form 

dF D 

-^-^(a;, Q 2 , xjp, t) = f(xjp, t) F 2 F ((3,Q 2 ) , (8) 

where as usual (3 = xjxjp and f(xp,t) is the pomeron "flux factor" which, if the diffrac- 
tively recoil system is a single proton, has the following form || [£| 

ft +^ AT B 2 {t) l-2a(t) / n \ 

f{xp,t) = N xjp , (9) 

where B(t) describes the pomeron coupling to a proton, see (|13|) f° r details, and iV is a 
normalisation factor. The function 



a(t) = a(0) + a' t (10) 

is the Regge trajectory of the pomeron. The model is based on the ("soft") pomeron 
exchange, which means that the intercept a(0) = 1.086 and the slope a' = 0.25GeV 2 . 
These values were obtained from the analysis of the high energy behaviour of total and 
elastic hadron-hadron cross sections (||10||). 

As we noted above F 2 F (/3, Q 2 ) is the pomeron structure function with the variable 
j3 playing the role of the Bjorken scaling variable. The factorization property (H) then 
follows as a direct consequence of the single pomeron-exchange mechanism of diffraction 
and the assumption that the pomeron is described by a single Regge pole. 
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The HERA measurements are found to satisfy factorisation property (||). Comparison 
of the observed values of the exponent n with that of formula after t integration, shows 
that the data confirms dominance of the "soft" pomeron exchange, since the "hard" or 
perturbative pomeron with intercept op(0) ~ 1.5 would give much steeper behaviour 
proportional to x~j^ or so. 

In the region of large Q 2 the pomeron structure function F<f(f3, Q 2 ) is expected to be 
described in terms of the QCD parton model which leads to a logarithmic violation of the 
Bjorken scaling, implied by perturbative QCD. In the leading logarithmic approximation 
the pomeron structure function is related in a conventional way to the quark qf{f3,Q 2 ) 
distributions in a pomeron: 

N f 

Ff(/3,Q 2 ) = 2 £e 2 /3gf(/3,Q 2 ) (11) 

i=l 

where are the quark charges and Nf is a number of active flavours (note that qf = qf). 
The quark densities qf(/3,Q 2 ) and the gluon density g F {(3,Q 2 ) of the pomeron evolve 
in Q 2 according to the Altarelli-Parisi evolution equations. However, it is necessary to 
provide input distributions at some reference scale Q 2 = Ql, since only the evolution 
is specified by perturbative QCD. The parton distributions in a pomeron at the initial 
scale Qq = 4GeV 2 were estimated in papers [|ll|]-[[li| using the "soft" pomeron interaction 
properties and " triple- Regge" phenomenology. In the following we present the results of 
the analysis of ref . [13| . 

First, it was assumed that the pomeron's flavour content is the same as that of the 
proton. Then, the quark distributions were parametrised in terms of a single function, the 
quark-singlet distribution, S JP (/5) = J2i P(lf(P) + Q~f((3))- The following parametrisation 
was found for Ql = AGeV 2 

E p {(3) = 0.0528 /r ' 08 (I -/3)+ 0.801 /3 (1 - /3) . (12) 

The "soft" term, proportional to /3 -0 08 , results from the triple- Regge limit of the diffrac- 
tive process ([!]) which applies when (3 — > 0, see Jl3| for details. The relatively small 



magnitude of the coefficient is a direct consequence of the small magnitude of the triple- 
pomeron coupling. The second term, which dominates for (3 ~ 1, results from an assump- 
tion that in this region of (3 the "soft" pomeron couples to a quark-antiquark pair with a 
form factor depending on a virtuality of the pair, see @ . 

The gluon distribution in a pomeron has a similar form 

p g F ((3) = 0.218 /r ' 08 (1 - (3) + 3.30 p (1 - (3) . (13) 

As before the "soft" term, dominating for small f3, results from the triple Regge limit. The 
analytical form of the second term, important for large (3, was assumed to be the same as 
for quarks and its coefficient was calculated assuming that the longitudinal momentum 
sum rule holds: 

jf dP (3 (£(<zf ((3) + qf(P)) + = 1 • (14) 
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The above parametrisation allows, via the evolution equations, the diffractive structure 
function (§) to be determined at any value of Q 2 accepted in perturbative QCD. The 
results are compared to HERA data, now given in terms of a new diffractive structure 
function F 2 D ^ defined as 

F 2 D(2 \P,Q 2 ) = r H F?®(/3,<?,x F )4x F , (15) 

Jx Pl 

where (3 is fixed during the integration. Therefore, assuming the simple factorisation 
formula (§), we see that the function F 2 D ^ is proportional to the pomeron structure 
function Ff. We find that generally the Bjorken scaling is observed, up to a possible 
mild (logarithmic) scaling violations. This fact favours a partonic interpretation of the 
pomeron structure function F 2 . 

In Figure 1 the solid curves show F 2 D ^ 2 \ obtained using the above parametrisation as 
the initial condition for Q 2 evolution. The comparison is done for HI and ZEUS data 
separately, because of different integration limits in flT5|) used by these experiments, see 
0. [l| for details. It should be noted that the increase of the computed pomeron structure 
functions F 2 D ^ with increasing Q 2 , up to the value (3 ~ 0.4, is related to the relatively 
large gluon distribution (|T3| ) in a pomeron which has a "hard" (1 — (3) spectrum at (3 ~ 1. 
This effect is nicely confirmed by the data, even for (3 = 0.65 the measured F 2 ^ rises 
with Q 2 . The above behaviour of the pomeron structure function F 2 ^ is in contrast to 
the behaviour of the proton structure F 2 (x, Q 2 ), which starts to decrease with increasing 
Q 2 already for x ~ 0.1. 

The presented parametrisation, however, cannot account for the rise of the diffractive 
structure function for the highest measured value of f3. The reason is simple, the relatively 
hard gluon spectrum is not hard enough. In order to understand this, let us consider the 
equation for the logarithmic slope of F 2 F , resulting from the Altarelli-Parisi equations 

where 5 = Y^fJi^f- The gluon term is always positive, whilst the first one could be 
negative, especialy for large (3 values, because of virtual emissions. Thus, F 2 ^ can rise 
with Q 2 for large (3 only when the gluon distribution is large enough at (3 ~ 1. 

The gluon distribution of the pomeron at large (3 is the most difficult density to 
estimate. In fact the analytical form of the second term in (|13|) is somewhat arbitrary 
and the momentum sum rule (|i~4]) used to determined its coefficient is also a source of 
controversy. The point is that the normalisation of the pomeron flux iV in ([8]) is arbitrary: 
any constant factor may be transferred from the pomeron flux to the pomeron structure 
function [irjj. Therefore, the magnitude of the momentum sum depends on convention. 



Lacking new theoretical ideas to estimate the gluon distribution of the pomeron, the 
only way to tackle the problem is to perform a QCD fit to the diffractive structure function 
data. 



4 



3 Pomeron parton distributions from QCD fits 



In order to clarify the gluonic content of a pomeron, we perform a global QCD fit to the 
diffractive data, published by the HI and ZEUS collaborations. A first attempt along this 



line was done in [15] using only HI data. We take the Regge factorisable form 

F^\x,Q\ XlP ) = x? Ff((3,Q 2 ) (17) 



where the partonic form QTT ) of is assumed. 

The Xp factor anticipates the dominant dependence of the pomeron flux ([|) on xp 
after the t integration. The exponent n and parameters of the parton densities at the 
starting scale Ql are allowed to vary simultaneously in the fit. Leading order Altarelli- 
Parisi evolution is used with four massless flavours. The parameter Aqcd in the running 
coupling constant a s is fixed in all fits. The momentum sum rule (|T3|) is not imposed, 
therefore, the normalization of the momentum sum, as well as the relative quark and 
gluon contributions, are determined by the fit itself. Only statistical errors are taken into 
account. 

The following forms of the singlet-quark and gluon distributions at the initial scale 
Ql = AGeV 2 are used 

S JP (/5) = A l {3 A > (1-(3) M , (18) 
/VG0) = B lf 3 B > (l-/3) B * , (19) 

whereas the nonsinglet quark distribution is set to zero. 

Table 1 summarises our results. In the first fit all seven parameters are allowed to 
vary, whilst Aqcd is fixed at 200MeV. The second fit shows how the parameters change 
if the value of n is fixed at the value determined by ZEUS. These two fits assume that all 
four flavours are massless, while in the last fit the charm flavour c is massive. In this case 

Ff(P> Q 2 ) = Ff {u ' d ' s \f3, Q 2 ) + Ff (c) ((3, Q 2 , m\) , (20) 

where F^ u,d,s ^ is given by QCD formula ([U]) with three massless flavours, and the charm 
contribution F 2 F ^ is driven by the gluon distribution function 



F 2 '(P,Q ,m c ) =e c / dy - C 2 -, — ) g (y,/i c ) , 21 

it Jap y \y Q 2 I 




where a = 1 + Am 2 c /Q 2 and /i 2 . = Am 2 c . The form of the Ci function as well as more 
details can be found in ref. |I2[ . The charm contribution is different from zero above 
the threshold for cc pair production, Q 2 (l — (3)/ (3 > 4m 2 . Not that there is no c quark 
distribution in this approach. 

The values of x 2 /dof confirms the observation made by HI and ZEUS that the present 
data are in good agreement with the factorisable form ([17]) of F 2 D . The fitted parameter 
n is in excellent agreement with the HI value (n = —1.19), although fixing it to the ZEUS 
value (n = —1.30) does not lead to a significant deterioration of the quality of the fit. 
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Parameters 


rll 1 


r lx z 


rll o 


— n 


1.18 


1.30 (fix) 


1.19 


A 

Ai 


0.066 


O.Ozl 


0.069 


A 


u.zy 


n 1 1 

U. 11 


n aa 


A 

^3 


n 79 

U. ( z 


n s/i 


n fin 


By 


1.22 


0.72 


1.16 


B 2 


3.13 


2.51 


5.00 


B 3 


0.31 


0.27 


0.000 


Aqcd 


200 


200 


200 


X 2 /dof 


114/96 


120/96 


110/96 



Table 1: The results of the fits to F 2 diffractive structure function data from HI and 
ZEUS experiments. Aqqd is fixed in all fits. The fitted parameters are described in the 
text. Only statistical errors are taken into account. 



Therefore, our analysis supports the idea of the dominance of the " soft" pomeron exchange 
in diffractive HERA data |fj, |4]] . The results of our studies are compared with the data in 
Figure 1. The dashed curves show from the massless flavour fit (Fit 1), whilst the 

dotted ones correspond to the fit with a massive c quark (Fit 3). The dotted curves at the 
bottom of each plot show the massive charm contribution Ff (c) to the if { 2) structure 
function in Fit 3 (the upper dotted curves). The curves from both fits give a reasonable 
description of the data, leading to a persistent rise of F 2 D(2) with log(Q 2 ), even for the 
largest value of (3. It is worthwhile to pointing out that the massive charm description 
offers an additional mechanism for the growth of F®^, especially for large (3 and for 
values of Q 2 well above the cc pair production threshold. 

Our singlet-quark and gluon distributions are plotted in Figure 2 as functions of j3 for 
different values of Q 2 . The continous curves correspond to the parametrisation presented 
in the previous section, the dashed and the dotted ones result from the massles and massive 
fits respectively. Their normalisation is determined by the choice iV = 2/tt in the pomeron 
flux, following the convention of refs. |6|, 0, |9|. With this normalisation the momentum 
sum rule (|14|) equals 1.7 for the parton distributions of the massles fit. The contribution 
to this sum from gluons is 85% at Q 2 = AGeV 2 , falling to 75% at Q 2 = 200GW 2 . 

The form of the gluon density at the initial scale Qq = 4 GeV 2 confirms the conclusion 
that the persistence of the rise of F 2 D(2) with log(Q 2 ) can only be reproduced by a large 
gluon (relative to quark) distribution at (3 ~ 1. However, it is worth repeating that the 
demand for such a gluon density cannot be demonstrated conclusively given the magnitude 
of the errors on the present data. In fact, within these errors, all our gluon distributions 
give a reasonable description of the HERA data. 
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4 Summary 



The analyses presented here show that the HERA diffractive deep inelastic scattering 
data can be described by a model in which the hard scattering of the virtual probe occurs 
on a colorless object with partonic structure. This object, a "soft" pomeron, appears 
naturally in the factorisable model of Ingelman and Schlein. The partonic distributions 
in the pomeron can be determined using the QCD parton model, where the Altarelli-Parisi 
evolution equations play an important role, especially for the determination of the gluonic 
content of a pomeron. Our analyses suggest that, within the present accuracy, a broad 
range of gluon distribution parametrisations is allowed. There is a strong indication, 
however, that the gluons constitute a large component of a pomeron, carrying most of a 
pomeron's momentum. It will be interesting to see whether more precise measurements 
confirm this indication. 

The potentially large gluon distribution can manifest itself in a significant heavy c 
flavour contribution to the diffractive structure functions and may also give rise to a large 
longitudinal structure function dF® /dxjpdt, which was neglected in the first analyses. 

More accurate measurements of the diffractive processes in deep inelastic scattering, 
over as wide a kinematical range as possible, will allow many important theoretically ques- 
tions to be addressed. The region of high (3 will clarify the gluon distribution problem. 
Measurements at very high y will allow the study of the longitudinal component of the 
diffractive cross section @, whereas measurements at xjp > 0.05 will test the pomeron fac- 
torisation picture, which may be violated due to lower lying reggeon exchanges (/, p, tt, ...). 
The universality of the extracted pomeron parton distributions may also be tested, by 
using them to describe exclusive hard diffractive scattering processes ]Tj| [T7|, |TS|| . 



In conclusion, deep inelastic diffraction is a promising subject for further experimental 
and theoretical studies in the coming years. 
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Figure 1: Comparison of the diffractive structure function F 2 D obtained in analyses of this 
paper with HI and ZEUS data. The solid curves correspond to the estimation described 
in section 2, while the dashed and dotted ones come from the massless (Fit 1) and the 
massive c quark (Fit 3) fits respectively. The dotted curves at the bottom show the charm 
contribution to the structure function obtained in the massive fit (the upper 

dotted curves). 
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Figure 2: The singlet and gluon distributions in a pomeron as functions of (3 for different 
values of Q 2 . The solid curves results from the estimation of chapter 2, the dashed 
and dotted ones are given by the massles (Fit 1) and massive (Fit 3) fits respectively. 
All distributions are normalised in accordance with the pomeron flux @ normalisation 
N = 2/tt. 
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